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Effect of Density Gradients in Confined Supersonic 
Shear Layers. Part II. 3-D Modes 

w 


s' 






5 


O , 

The effect of basic flow density gradients on the supersonic wall modes were 
investigated in Part I of this analysis. In that investigation only the 2-D modes 


were studied. Tam and Hu investigated the 3-D modes in a confined vortex 


sheet and reported that the first 2-D Class A mode (Aqi) had the highest 


growth rate compared to all other 2-D and 3-D modes present in the vortex 

-4- 

sheet for that particular set of flow parameters. They also showed that this 7* jE 
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result also held true for finite thickness shear layers with < 0.125. For free 


shear layers, Sandham and Reynolds showed that the 3-D K-H mode became 


the dominant mode for Me > 0.6. Jackson and Grosch investigated the effect ►- 


of crossflow and obliqueness on the slow and fast modes present in a M c > 1 
environment and showed that for certain combination of crossflow and wave 
angles the growth rates could be increased by up a factor of- 2 with respect 
to the 2-D case. Two major differences between the confined and unconfined 
cases are the type of modes present in the two cases and the restriction on the 
spanwise wavenumber in the confined case. The modes present in a confined 
supersonic 2-D shear layer are of the acoustic type (wall modes) as opposed to 
the vorticity modes present in the free shear layer, and these wall modes have 
higher growth rates than the vorticity modes. Also, due to the side walls, the 
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confined shear layer only has discrete spanwise wavenumber solutions. 

The case studied here is a confined shear layer shown in Part I. All solution 
procedures and basic flow profiles are the same as in Part I. The effect of density 
gradients on the 3-D modes present in the density ratios considered in Part I 
are investigated. 

II. Flow Model and Governing Equations 

As in Part I, the model which is studied here is a confined compressible 
shear layer formed by two gases with different velocities, densities and prop- 
erties, but with the same constant pressure. Figure 1 shows the configuration 
used in this analysis. The subscript 1 is used for the quantities related to the 
high speed freestream and the subscript 2 for the quantities of the low speed 
freestream. The streamwise coordinate is x, the spanwise coordinate is 2 , and 
the cross-stream coordinate is y. The aspect ratio of the channel is taken to 
be 2 ( B * = 2 H*). For simplicity, free slip wall boundary conditions are as- 
sumed at the walls of the channel. It is also assumed that the flow is inviscid, 
non-conducting, and non-diffusive. For this situation the governing equations 
are the Euler equations for a two species system. 

The equations are non-dimensionalized using the fast (upper) freestream 
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quantities, p*, U*, T* and the height of the channel, H*. Thermodynamic 
properties are also non-dimensionalized by the upper freestream thermody- 
namic properties, Cp Q and . Based on this non-dim ensionalization, the den- 
sity ratio ^ and the velocity ratio are defined as p 2 and U 2 , respectively. 
Once the equations are non-dimensionalized, they are linearized around a par- 
allel basic flow ( p(y ), P , U (y), C a (y ))■ These basic flow quantities can be found 
by solving the compressible boundary layer equations using a similarity vari- 
able for Pr = Sc = Le == 1. This procedure was discussed in Part I. 

Normal modes are assumed for these infinitesimal disturbances with the 
form: 


q = q(y) exp[i(kx + fiz — wt)] (2.1) 

where q is the eigenfunction, k and /3 are the streamwise and spanwise wave 
numbers and u> is the frequency. In general k, (3 and u> are complex. Once these 
normal modes are substituted into the equations of motion, a single O.D.E. 
can be found for the disturbance pressure eigenfunction: 
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where 



_ _ C Po Cg + C p|> ( 1 ~ Cg) 

C va Ca “I" c *t(l — Ca) 

and Mi and 71 are, respectively, the Mach number and the ratio of specific 
heats of the high speed freestream. 

The above equation is solved subject to the boundary conditions 

Dp = 0 y = (2.3a) 


Dp = 0 z — ±1 (2.36) 

which come from the y and z-momentum equations by setting the normal ve- 
locity at the walls to zero. It is the latter boundary condition that leads to a 
discrete spectrum of spanwise wavenumbers. Therefore the spanwise wavenum- 
ber becomes: 


/ 3 — ran 


(2.4) 
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for an aspect ratio of two, and the disturbance pressure takes the form: 


p (x, y , z, t ) = p(y ) cos(mnz ) exp[i(kx — wf)] (2.5) 

t 

Outside the shear layer the vortex sheet solutions hold. They can be written 
as: 


p(y) = Acos[Ai(l — y)] 6<y<l (2.6 a) 

p(y) = 5cos[A2(1 + y)] - 6 > y > -1 (2.6 b) 

where 

Aj = [M?(u> - kUi ) 2 - k 2 - (m7r) 2 p (2.7a) 

A 2 = [M 2 (4)V - kU 2 ) 2 - k 2 - (mn) 2 ]> (2.7 b) 


Following Mack, we define a 3-D wave number as: 


k = (fc 2 + (m7r) 2 )^ (2.8) 


one can also define a complex convective Mach number, M c (y) as: 
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( 2 . 9 ) 


M c (y) = 


rpyi iMf 

7& 2 


(w - fct /) 2 


therefore equation (2.2) can be written as: 




( 2 . 10 ) 


This form of the O.D.E. for the disturbance pressure eigenfunction clearly 
shows the role of the parameter M c and its significance in the type of solutions 
possible. Now, outside the shear layer for a neutral wave this convective Mach 
number is a constant real value, i.e. (M c i = M c i, M c 2 = M c 2 ), 


M c 1 = Mi 


kU - 


u> 


= M\cos^(l - C p hx) 


(2.11a) 


and 


M c 2 = Mi cosi){U2 ~ Cphx) (2.116) 

where Cphx = * and 0 is the wave Angle. As shown for the 2-D case, 

when M c2 1 non-inflectional neutral solutions of the type Cphx — 1 
possible and when M c \ > 1 then non-inflectional neutral solutions of the type 

Cphx = ^2 are possible. 
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Based on the above definitions and the results obtained in Part I, the 


following cases are investigated: 


Case/ : p 2 — 1.398 m — 1 

Casell : p 2 = 1-398 m = 3 

Caselll : p 2 = 3.000 m = 1 

CaselV : p 2 = 3.000 m = 3 

In all the cases, the fast stream gas is He and the slow stream gas is JV 2 . 
Also, the convective Mach number is 1.836, and the velocity ratio is 0.276. 

Results and discussion 

Case I 

As shown in Part I, the basic flow profiles have three generalized inflection 
points giving rise to three separate types of modes. These modes are all 
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supersonic wall modes and the necessary condition for their existence is a 
trapped region of supersonic flow with respect to their phase speeds. In order 
to use this condition, the maximum convective Mach number in the channel, 
\Mc(y)\max , is plotted in the frequency and real wavenumber plane. Figure 2 
shows the contours of this maximum convective Mach number. The dashed 
lines correspond to \Mc(y)\max < 1. Any modes which are present within 
the ‘wedge’ ( \M c (y)\max = 1) are subsonic with respect to both streams and 
therefore are not acoustic in nature. Figure 3 shows the real wavenumber vs. 
frequency of all the modes present when the spanwise wavenumber is fixed to 
= 7 r. The fine dashed line in the form of a wedge is the \M c (y)\max = 1 
line. Two modes originate in that wedge and are subsonic with respect to 
both streams until they exit the wedge. Figure 4 shows the phase speed of the 
modes present for this case. Several differences arise when we compare Figure 
4 to Figure 6 in Part I. First, we now have a mode present in the limit as 
the frequency goes to zero. This mode is originally subsonic, and as it turns 
supersonic with respect to the fast stream it is very much like the combination 
of the first three 2-D Class B modes (upper portion of the phase speed curves 
in Figure 6 in Part I). This mode will be called the D\\ mode. The first Class 
C mode C\\ also starts from within the wedge described above, but its phase 
speed exits the wedge very quickly and only slightly changes the characteristics 
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of the mode. The rest of the Class A, B and C modes are all entirely acoustic 
modes and their phase speeds are similar to those in the 2-D case. One must 
note that the first two Class A and C modes are different in the sense that 
each has ‘picked up’ the other’s phase speed curve. Figures 5 a,b,c and d show 
the growth rate of the Class A,B,C and D modes respectively. The growth 
rate curve for the A \ 2 mode has only one peak and its maximum growth rate 
is a few percent higher than than that of the 2 -D mode. The growth rate curve 
for the A 13 mode is almost identical to its 2-D counterpart despite the fact 
that the wave angle at the maximum growth rate of for this mode is about 30 
degrees. Mode A\\ has very small growth rates (as do modes Bn, B 12 , B 13 ) 
and is not shown. Since the Class A and C modes have switched portions of 
their phase speed curves with one-another, one must also discuss the growth 
rates of the Class C modes with conjunction to the Class A modes. Instead 
of identifying the 2 -D and 3 -D counterparts of each growth rate peak, one can 
make a general statement that except for the A i2 , in the region where the 
waves are 3 -D, the growth rates are lower than those found for the 2 -D case. 
In fact mode C\\ has a growth rate which is about 25% lower than its 2 -D 
counterpart. The growth rate curves of the class B modes are very much like 
the 2 -D modes. The D\\ mode present for this spanwise wave number has 
one subsonic growth rate peak and three supersonic peaks. The supersonic 
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portion of the phase speed curve for this mode is very much like the curve for 
the I?oi mode for a density ratio of about 1.65 where the first three Class B 
modes have undergone a resonant interaction. In general, if one plots all the 
growth rates of the 2-D and 3-D modes (/3 — 7r), ope will find that the Cq\ has 
the highest growth rate. 

Case II 

In this case, the density ratio is the same as above, but the spanwise 
wavenumber is now equal to 37 t. Figure 6 shows the maximum convective 
Mach number contours. The wedge corresponding to the \M c (y)\max 1 
contour is much larger than that in Case I. This implies that the region over 
which purely subsonic modes can exists is much greater. Figure 7 shows the 
real wavenumbers of the modes present when the spanwise wave number is 
set to 37r. Now both modes which start within the subsonic wedge have zero 
frequency limits and both will be labeled as Class D modes. Figure 8 shows 
the phase speed of the modes for this case. There are Class A,B and C modes 
with very small growth rates within this frequency region and are not shown. 
Figure 9a shows the growth rates for some of the Class A and B modes and 
Figure 9b for the Class D modes. Very much like Case I, it can be concluded 
that the 3-D modes have lower growth rates than the 2-D modes. Therefore 
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for the test case parameters the 2-D, the Cq\ mode has the highest growth 
rate. Tam and Hu also found similar results for their v4qi mode. 

Case III 

This case corresponds to a density ratio of 3.0 for which the basic flow 
profiles only have one generalized inflection point. Thus, as shown in Part I, 
only Class A and B modes are present in the 2-D analysis. Figure 10 shows 
the real wavenumber for the modes present. Only one mode exist in the 
subsonic wedge and is present in the zero frequency limit. The phase speed 
vs. frequency is given in Figure 11 and Figures 12a and b show the growth 
rates of the Class A modes and the D\\ mode respectively. The 3-D Class 
A modes for this case have smaller growth rates than their 2-D counterparts. 
Very much like Case I, the D\\ modes’ supersonic continuation is very similar 
to the £ 0 i mode present for these parameters . Figure 13 shows the growth 
rates of both modes on the same graph. It can be seen from this figure that 
the maximum growth rate peak for the 3-D mode is slightly higher than the 
2-D mode. More about the growth rate of this mode will be said in the next 
section. 

Note that the Class B modes were not shown due to fact that they possess 
very small growth rates. 
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Case IV 


The parameters for this case are the same as in Case III, however the span- 
wise wavenumber is increased to 37T. Figures 14 and 15 show the wavenumber 
and phase speed of the modes present. Again, like Case II the sonic wedge is 
much greater than for the smaller wavenumber (Case III). Figure 16 and 17 
show the growth rates of the Class A and the D\\ modes. The Class A modes 
have much lower growth rates than their 2-D counterparts. However, the D\\ 
mode a higher growth rate (by 5 %) than the .Boi mode. Also, the frequency 
for the maximum growth rate has shifted to the left. In order to investigate 
this effect further, the mode D\\ for (3 = is calculated and its growth rate 
along with the D\\ modes for (3 — n and 37r and the Bqi modes are plotted in 
Figure 18. It is clear from this figure that the rise in the growth rate due to the 
larger growth rates in the subsonic regions of these modes eventually stops and 
actually the growth rates start to decrease beyond a spanwise wavenumber of 
57T. 

Based on all the results shown for all the cases considered, one can conclude 
that 3-D effects are not that important in the sense that they don’t give rise 
to much larger growth rates as is the case in free shear layer. However, they 
do give rise to modes which have almost the same phase speeds and growth 
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rates as the 2-D modes which can give rise to non-linear resonant interaction 
in the evolution of the full 3-D shear layer. 
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Growth Rate vs. Frequency for 3-D Class B Modes 
(B = 2H, p = 2mn/B, m = 1) 
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Phase Speed vs. Frequency for All 3-D Modes 
Based on the Boundary Layer Profiles 
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Growth Rate vs. Frequency for the D 01 Mode 
(B = 2H, (3 = 2 rri 7 c/B, m = 3) 
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